﻿Nicolae Cotfas ZS INTRODUCERE IN ANALIZA FOURIER-LAPLACE - Versiunea 10 martie 2015 - http://fpcm5 fizica unibuc ro/~ncotfas/ E-mail: ncotfas@yahoo com Tel: 074 278 4634 Introducere Prezentare centrată pe exemple, aplicații și calcul numeric cu teoria redusă la elementele strict necesare, prezentate accesibil, precis, clar și concis București, 2013 Nicolae Cotfas 3 Cuprins 1 Introducere 7 1 1 Elemente de trigonometrie 7 1 2 Numere complexe 10 1 3 Funcții complexe de variabilă reală 13 1 4 Derivarea funcțiilor complexe de variabilă reală 16 1 5 Integrarea funcțiilor complexe 19 1 6 Serii de puteri 23 1 7 Spații cu produs scalar 25 2 Serii Fourier 31 2 1 Exemple 31 2 2 Aplicații 50 2 3 Definitii si rezultate fundamentale 50 3 Transformarea Fourier discretă 55 3 1 Definitie si exemple 55 3 2 Transformarea Fourier rapida 63 3 3 Proprietati ale transformarii Fourier 64 3 4 Functii proprii ale transformării Fourier 64 3 5 Transformarea Fourier bidimensionala 65 3 6 Aplicatii 70 3 7 Definitii si rezultate fundamentale 70 4 Transformarea Fourier a funcțiilor 77 4 1 Exemple 77 5 4 2 Aplicații 82 4 3 Definiții și rezultate fundamentale 82 5 Transformarea Fourier a distribuțiilor 89 5 1 Exemple 89 5 2 Aplicații 101 5 3 Definiții și rezultate fundamentale 101 6 Transformarea Laplace a funcțiilor 119 6 1 Exemple 119 6 2 Aplicatii 119 6 3 Definitii și rezultate fundamentale 119 7 Transformarea Laplace a distribuțiilor 121 7 1 Exemple 121 7 2 Aplicatii 121 7 3 Definitii și rezultate fundamentale 121 8 Transformarea Laplace finită 123 8 1 Exemple 123 8 2 Aplicatii 123 8 3 Definitii și rezultate fundamentale 123 6 Capitolul 1 Introducere Pe parcursul acestui capitol vom prezenta cateva elemente ale apartului matematic care va fi utilizat în capitolele urmatoare 1 1 Elemente de trigonometrie 1 1 1 Se stie ca unghiurile posibile, masurate în radiani, pot fi reprezentate in mod natural pe circumferința cercului trigonometric (vezi fig 1 1) Pentru fiecare t E R, unghiurile t — 4n, t — 2n, t, t+2n, t+4n, se reprezinta în acelasi punct 1 1 2 Definiție Functia cosinus este functia (vezi fig 1 1) cos : R —> [—1,1] : t cos t unde cos t este coordonata proiectiei pe axa orizontala a punctului în care se reprezinta t pe circumferinta cercului trigonometric 1 1 3 Definiție Functia sinus este functia (vezi fig 1 1) sin : R —> [—1,1] : t sin t unde sin t este coordonata proiectiei pe axa verticala a punctului în care se reprezinta t pe circumferinta cercului trigonometric 7 8 Aplicații ale transformărilor Fourier și Laplace Figura 1 1: Funcțiile cosinus și sinus 1 0 0 5 5 -10 -5 10 -0 5 Figura 1 2: Graficul funcției cosinus 1 1 4 Direct, din definițiile anterioare rezultă că: - Funcțiile cosinus și sinus sunt periodice cu perioada principală 2-tf cos(t + 2-tf) = cost sin(t + 2-tf) = sint oricare ar fi tGR - Funcția cosinus este funcție pară și sinus impară cos(—t) = cost sin(— t) = — sint oricare ar fi t G R - Funcțiile cosinus și sinus verifică relația fundamentală cos21 + sin21 = 1 oricare ar fi tGR 1 1 5 Unele valori se pot determina simplu utilizând geometria elementară t = 0 t — - b “ 6 t — - t — - t — -b “ 2 t = 7F j- 3tf b ~ 2 sint 0 1 2 V2 2 V3 2 1 0 -1 cost 1 ^3 2 V2 2 1 2 0 -1 0 (1-1) Introducere 9 -5 1 0 0 5 10 -0 5 -1 0 Figura 1 3: Graficul funcției sinus 1 1 6 Utilizând geometria elementară se poate arăta că au loc relațiile fundamentale sin(a ± b) = sin a cos b ± cos a sin b oricare ar fi a, bGR cos(a ± b) = cos a cos b =f sin a sin b adică sin(a + b) = sin a cos b + cos a sin b sin(a — b) = sin a cos b — cos a sin b oricare ar fi a, 6eR (1 2) cos(a + b) = cos a cos b — sin a sin b cos(a — b)= cos a cos b + sin a sin b 1 1 7 Utilizând (1 1) și (1 2) se pot obține relațiile șina = cos (^ — a) cos a = sin (^ — a) sin 2a = 2 sin a cos a cos 2a = cos2 a — sin2 a = 2 cos2 a — 1 = 1 — 2 sin2 a cos2 a = 1+c°s2a sin2 a = x~c°s2a șina cosb = |[sin(a+&) +sin(a—&)], cosa cosb = |[cos(a+&) +cos(a—b)], șina sinb = |[cos(a—b) — cos(a+&)] sin a + sin b = 2 sin cos sin a — sin b = 2 sin cos cos a + cos b = 2 cos cos cos a — cos b = — 2 sin sin 10 Aplicații ale transformărilor Fourier și Laplace 1 2 Numere complexe 1 2 1 Definiție In cazul unui număr complex z = x + yi Re z = x se numește partea reală a lui z Im z = y se numește partea imaginară a lui z z = x — yi se numește conjugatul lui z |z| = \/x2 + y2 se numește modulul lui z 1 2 2 Direct din definiție rezultă următoarele relații zi ± z2 = zi ± z2 |z| = |z| Re z = — ziz2 = zi z2 |z|2 = zz Im z = (zn) = (z)n (z) = z z = Re z+i Im z 1 2 3 MATHEMATICA Re[x+y I], Im[x+y I], Abs[x+y I], Conjugate[x+y I] In :=I In :=Sqrt [-4] In : = (3+2 I)*2 In : = (3+2 I)/(5-I) Out = ii Out =2i Out =5+12i Out = 2 + 2 In :=Re In :=Im In :=Abs In :=Conjugate Out =3 Out =4 Out =5 Out =3-4i 1 2 4 Se stie că i2 = —1 si (a + i) ± (a + ^i) = (a ± a) + (b ± ^)i (a + bi) (a + ^i) = (aa — b^) + (a^ + ba)i a | bi (a+bi)(a-Și) (aa+bȘ)+(-aȘ+ba)i a+Și (a+Și)(a-Și) a2+Ș2 ’ De exemplu, (2 + 3i) + (1 + 2i) = 3 + 5i (2 + 3i) — (1 + 2i) = 1 + i (2 + 3i)(1 +2i) = —4 + 7i 2|3i (2|3i)(1-2i) 8-i 8 1i 1+2i (1|2i)(1-2i) 5 5 5i’ Introducere 11 1 2 5 Teoremă Relațiile a) |(xi+yii) + (X2+ysi)| x/x2 = |x| |x + yi| \ x2 + y2 > yfy2 = |y| iar relatia 12 Aplicații ale transformărilor Fourier și Laplace \ x2 + y2 0 converge la a+bi si scriem dacăa lim (xn+yni) = a+bi n^^ lim |(xn + yni) - (a + bi)| = 0 n^^ 1 2 9 Teoremă Avem lim (xn+yni) = a+bi C, f (t) = u(t) + i v(t) definite pe un interval I C R functia u : I —> R 14 Aplicații ale transformărilor Fourier și Laplace reprezintă partea reală a lui f, iar v : I —> R partea imaginară a lui f 1 3 2 Relația f : R —> C, f (t) = e14 = cos t + i sin t definește o functie periodica f (t + 2n) = cos(t + 2n) + i sin(t + 2n) = cos t + i sin t = f (t) cu partea reala u : R —> R, u(t) = cos t Si partea imaginara v : R —> R, 1 3 3 Din relatiile e14 cos t + i sin t v(t) si = sin t e 14 = cos t — i sin t se obtin formulele e14 + e-14 cos t = 2 si e14 e—14 sin t = 2i 1 3 4 Pentru T E (0, rc>) fixat, functia complexa f (t) = e1'4 f : R —> C, cos t + i sin 2^ t este o functie periodica cu perioada principala T , ; 2n ii | rri\ ; 2n 4 ; 2n 4 , f (t + T) — e1 t (4+T) — e1 t 4 e12n — e1 t 4 — f (t) f (t 1 J) — e ' — e ' e — e ' — f (t)- 1 3 5 Teoremă Fie t° E I și f : I —> C, f (t) = u(t)+i v(t) Avem {u este continuă în t° v este continuă în t° Demonstrație Din relatia u(t°\l r f(t1-)> f(t1+)> ••• , f(tn-), f(tn+), f(b—) exista si sunt finite Spunem ca funcția f ca este continua pe porțiuni în R dacă este continua pe portiuni în orice interval [a, b] 1 3 8 Functia (vezi Fig 1 5) f : [-1, 2] -^ R, t2 - 1 dacăa dacăa -1 R este o funcție deriva,bilă in punctul toEl dacă există și este finită limita numită derivata, lui f in t0 1 4 2 Exemple a) Funcția / : R —> R, f(t) = t3 este derivabilă in orice punct to GR deoarece lim j|m £ -£) j|m^2 +t'o) = 3 t'o- t->-a t — to x^a t — to In acest caz, fțto) = 3 tg, oricare ar fi toGR, adică avem (t3)' = 3t2 b) Funcția /: [0, oo) —>R, f(t) = x/i este derivabilă in orice punct to € (0, oo) iim LfLzAA t — tg t — tg = lim —= — x/t + x/to 1 In acest caz, /'(to) = oricare ar fi to G (0, oc), adică avem (Vt)' = 1 4 3 Derivatele unor funcții uzuale Introducere 17 Funcția Derivata Domeniul Condiții f R —— R f (x) = c f '(x) = 0 R f R —— R f (x) = xn f z(x) = nxn—1 R n G N* f (0, w) —— R f (x) = xa f z(x) = axa—1 (0, w) a G R f R* —— R f (x) = 1 f/(x) = — x R* f [0, w) —— R f (x) = y/x f/(x) = 2-x (0, w) f [0, w) —— R f (x) = n/x f '(x> = n -1 (0, w) n G 2N* f R —— R f (x) = /x f ' = n 'Jxn—i R* nG2N+1 f (0, w) —— R f (x) = ln x f/(x) =1 (0, w) f R —— R f (x) = ax f z(x) = ax ln a R 0 C, f (t) = u(t) + i v(t) derivabile în to avem lim f(t) - f ' = jim u(t) - ' + i lim v(t) — t^to t — to t — to t — to adica, f '(to) = u'(to) +i v' (to) 1 4 7 Pentru a GR fixat, avem ie1'" / = (cos at) + i(sin at) = —a sin at + i a cos at = i a(cos at + i sin at) = ia eiat adicăa, = ia eiat 1 4 8 Daca :R —>R este o funcție derivabilă atunci (e^(t)) = (cos ^>(t))' + i (sin ^>(t))' = — ^>'(t) sin tp(t) + i ^>'(t) cos ^>(t) = i ^'(t)(cos R sunt doua primitive ale funcției f: IR atunci exista c E R astfel încat Fi = F2 + c Vom nota cu C multimea functiilor constante definite pe intervalul considerat Multimea primitivelor unei funclji f se notează cu f f (x)dx, adica y f (x)dx = { F: IR | F este primitivă a lui f } 1 5 3 Primitivele unor funcții uzuale (I este un interval inclus în domeniul maxim de derivabilitate al primitivelor) Funcția Mulțtimea primitivelor Intervalul Condiții f (x) = 1 f dx = x+C IC R f (x) = xn f xndx = xn—1 + C IC R n E N f (x) = xa J x"'dx = o+î xa+1 + C IC (0, ) a E R — {—1} f (x) = 1 J 1 dx = ln |x|+ C IC R—{0} f (x) = ex J ex dx = ex +C I c R f (x) = ax J ax dx = ax +C I c R 0 0 f (x) = , 2 2 J \ J v/x5+a5 f , 5 5 dx = ln(x+Vx2+a2 ) + C J vx5+a5 \ / IC R a=0 f (x) = a2 + x5 f -dx = 1 arctgx + C IC R a=0 f (x) = x5-a5 f -21^ dx = 21 ln x +C d x5-a5 2a x—a IC R—{±a} a=0 f (x) =sh x f shxdx = chx+C IC R f (x) =ch x ch x dx = sh x+C IC R 20 Aplicații ale transformărilor Fourier și Laplace 1 5 4 MATHEMATICA: Integrate[f[x], x] In :=Integrate[f[x], x] i—— Out = J f (x) dx In :=Integrate[x“a, x] i—— Out = xj+- In :=Integrate[a“x, x] i—— Out = L In :=Integrate , x] i—— Out =Log [x+Va2 +x2 ] 1 5 5 Teoremă (Integrarea prin părți) Dacă funcțiile f,g: I—R sunt derivabile și au derivate continue atunci y f (x) g'(x) dx = f (x) g(x)- f '(x) g(x) dx 1 5 6 Exercițiu Să se calculeze y"xex dx, y ln xdx Rezolvare Avem J xexdx = f x(ex)'dx = xex — f exdx = xex — ex + C f ln x dx = f x' ln x dx = x ln x — f x (ln x)'dx = x ln x — f dx = x ln x — x + C 1 5 7 MATHEMATICA: Integrate[f[x], x] In :=Integrate[x Exp[x], x] i—— Out =ex(- 1+x) In : =Integrate [Log [x] , x] i—— Out = -x+x Log [x] 1 5 8 Teoremă (Schimbarea de variabila) Fie I, JC R intervale și I —— J ——— R două functii Dacă ip : I —— J este derivabilă și F este o primitivă a lui f atunci I—R:x — (Fo^)(x) = F(^(x)) este o primitiva a functiei I—R: x — f (^(x)) p'(x), adica dt) o^> Demonstratie Avem dxF(^(x)) = F'fp(x)) p'(x) = f(^(x)) p'(x) y f (^(x)) ^'(x) dx = U f (t) 1 5 9 Exercițiu Să se calculeze Rezolvare Avem dx = 2 f 2ț^țdx = 2 f e^^Vx/dx = 2 (f e*dt)t== 2^® + C Introducere 21 1 5 10 MATHEMATICA: Integrate[f[x], x] In :=Integrate[Exp[Sqrt[x]]/Sqrt[x], x] i—— Out =2ev'x 1 5 11 Funcția complexă F : I —— C, F(t) = U(t) + i V(t) este o primitiva a funcției f : I —— C, f (t) = u(t) + i v(t) daca Fz(t) = f (t) oricare ar fi t GI adicaă dacăa Uz(t) + i Vz(t) = u(t) + i v(t) oricare ar fi t GI Ultima relație are loc daca si numai daca U'(t) = u(t) /ATI C* Q TO or TI oricare ar V'(t) = v(t) t GI adica daca U este primitiva a lui u si V este primitiva a lui v 1 5 12 Pentru a = 0 fixat, F : I —— C, F(t) = — eiat ia este o primitivăa a functiei f : I -— C, f (t) = eiat deoarece / i V I -1- e^t ) eiat ia 1 5 13 Știm ca dacă F:I —— R este o primitiva a functiei f: I —— R atunci f f (t) dt = F(t)|t=O = F(0) - F(a) J a oricare ar fi intervalul [a, 0] CI 1 5 14 Dacă F : I —— C, F(t) = U(t) +i V(t) este o primitiva a functiei f : I —— C, f (t) = u(t) + i v(t) atunci ,i' ' f (t) dt = /O u(t) dt + i /O v(t) dt = U(t)\t=a + i V(t)\t=a = F(0) - F(a) oricare ar fi intervalul [a, 0] CI 22 Aplicații ale transformărilor Fourier și Laplace 1 5 15 Exercițiu Să se arate că eint elkt dt = Snk Rezolvare Daca n = k atunci 1 r n eint eint 2n J-n •n dt = 1 •n iar dacă n = k atunci — / eint eikt dt = 2n J-n — r ei(k-n^t dt= 2n J-n 2ni(k — n) ei(fc-n)t n = 0 1 —n 1 5 16 O funcție continua pe porțiuni este integrabila pe orice interval finit 1 5 17 Teoremă (Formula de integrare prin parti) Daca funcțiile f,g: I—R au derivate continue pe intervalul I atunci I f,(t) g(t) dt = f (t) g(t)\ba —[ f(t) g,(t) dt Ja Ja oricare ar fi a,b EI Demonstrație Utilizand formula Leibniz-Newton obtinem f (t) g(t)ia = fOb(f ' g)'(t) dt = fb(f '(t)g(t) + f (t)g'(t))dt aa a = Jaf/(t) g(t) dt + tâf (t) g/(t) dt- 1 5 18 Exercițiu Sa se calculeze integralele yn yl / x cosxdx / x2ex dx Jo Jo Rezolvare Utilizand integrarea prin parti obținem fo t cos tdt = /fit (sin t) dt = x (sin t)|J — ffisin tdt = cos tio = —1 — 1 = —2 f 1t2 et dt = f 1t2 (etj) dt = t2 et|1 — 2 f 1t et dt = e — 2 f 1t (et)/ dt 0 e = 0 (e ) = e 0 — 0 e = e — 0 (e ) = e — 2t et|1 + 2 L1et dt = —e + 2 et|1 = e — 2 I u •j u IU 1 5 19 MATHEMATICA: Integrate[f[t], {t, a, b}] NIntegrate[f[t], {t, a, b}] In :=Integrate[t Cos[t], {t, 0, Pi}] i—> Out =-2 In :=Integrate[t“2 Exp[t], {t, 0, 1}] i—— Out =-2+e Introducere 23 1 5 20 Teoremă (Schimbarea de variabilă) Fie funcțiile [a, b] —J -—-+ R, unde J C R este un interval Daca f este continua și f este derivabila cu derivata continua atunci rb , r FF) / f (F(t)) F (t) dt = f (x) dx Ja J R este functie para, adica f (—t) = f (t) oricare ar fi t, atunci J“ f (t) dt = J0a f (t) dt + Joa f (t) dt = f '' f (—t) dt + J0a f (t) dt —a —a 0 —a 0 = — Ja0 f (t) dt + J01 f (t) dt = J01 f (t) dt + J01 f (t) dt = 2 Ja f (t) dt- 1 6 Serii de puteri 1 6 1 Stim ca nu exista dacaă x V 1 xn+1 = lim 1—x n^tt 1 1—x 1—x adică funcția (-1,1) —> R : 1 x - 1-x admite dezvoltarea in serie de puteri 1 1- — = 1 + x + x2 + - x 1 6 3 In cazul unei serii de forma (numita serie de puteri) a0 + a1(x — x0) + a2(x — x0)2 + a3(x — x0)3 + • • • exista 0 R astfel încat a) (x,ay+fiz) — a(x,y) + 0(x,z), ^x,y,zG V si Va,^ GR b) (x,y) — (y,x), Vx,y G V c) (x, x) > 0, Vx G V si (x, x) — 0 x — 0 1 7 2 Exemplu Pe spatiul vectorial real R2 — {x — (xi,x2) | xi,x2 G R } unde (xi,x2) + (yi,y2) — (xi +yi,x2 + y2) a(x1,x2) — (axi,ax2) relatia (x,y) — xiyi + x2y2 defineste un produs scalar Numărul ||x|| — \/ (x,x) adicăa ||(xi,x2)|| — y/ x2 + x2 reprezinta lungimea vectorului x — (xi,x2) 26 Aplicații ale transformărilor Fourier și Laplace Figura 1 6: Lungimea vectorului x = (x1,x2) din R2 1 7 3 Utilizand formula cos(a — b) = cos a cos b + sin a sin b obținem relația (a se vedea figura 1 7) (x, y) = xiyi + x2y2 = ||x|| ||y|| (cos a cos b — sin a sin b) = ||x|| ||y|| cos(a — b) din care rezulta ca produsul scalar (x, y) a doi vectori x si y este egal cu produsul lungimilor vectorilor înmultit cu cosinusul unghiului dintre ei In particular, vectorii sunt ortogonali (perpendiculari) daca si numai daca (x,y) = 0 Figura 1 7: Produsul scalar a doi vectori din R2 1 7 4 Definiție Prin produs scalar pe un spatiu vectorial complex V se întelege o aplicatie Introducere 27 (,} : V x V —> C astfel încât a) (x,ay+ftz) = a(x,y) + fl(x,z), Vx,y,zGV si Va,ȘGR b) (x,y) = (y,x), Vx,y G V c) (x, x) > 0, Vx G V si (x, x) = 0 x = 0 1 7 5 Exemplu Pe spațiul vectorial complex C2 = {x = (xi, x2) | xi,x2 G C } unde (xi,x2) + (yi,y2) = (xi +yi,x2 + y2) a(x1,x2) = (axi,ax2) relatia (x,y) = xi yi + x2 y2 defineste un produs scalar Numărul ||x|| = \/ C | este funcție continua } unde W + W(t) = W + WL W)(t) = a W) relatia 1 fn (M) = -/ (x) = 2 a0 + a1 cos t + b1 sin t ao = 1 f-n +(t) dt ai = 1 +(t) costdt bi = n +(t) sin tdt și în general, 1 k +(x) = - a0 + (an cos nt+bn sin nt) 2 n=1 an = 1 +(t) cos nt dt bn = n J-n +(t) sin nt dt Serii Fourier 33 Rezolvare Avem n f-n ^(t) dt = n f-n a0 + ai cos t + b1 sin dt = a0 1 ^(t) cos t dt = 1 ^2a0 + a1 cos t + b1 sin t^ cos t dt = a1 1 ^(t) sin t dt = 1 (2a0 + a1 cos t + b1 sin t^ sin t dt = b1 2 1 6 Dacă seria este convergentă si poate fi derivată termen cu termen atunci 1 ~ ( an = 1 f-n C cu perioada 2n pentru care coeficienții 1 [n 1 [n an = — ^>(t) cosntdt, bn = — ^>(t) sinntdt n J—~ n J—~ există i se asociază seria 1 00 -a0 + ^2(an cos nt+bn sinnt) 2 n=1 (2 2) numita seria Fourier (trigonometrică) a lui 2 1 8 Daca notăm funcțiile din sirul (2 1) cu e0, e1, e2, atunci (en, 6k) = dnk coeficiențtii seriei Fourier verificaă relațtiile a0 = (60 , C periodica cu perioada 2n Pentru ca o functie sa coincida cu suma seriei Fourier asociate trebuie ca ea sa fie periodica cu perioada 2n 2 1 13 Functiile de forma : [—n,n) —> C pot fi identificate cu funcțiile periodice : R —> C cu perioada 2n (obțnute folosind prelungirea prin periodicitate) 2 1 14 Fie este o funcție periodică cu perioada 2n Oricare ar fi t0 ER avem n f-n ^(t) cos ktdt = n ^(t) cos ktdt+n c cos kt dt = n J-°7T ^(t + 2n) cos k(t + 2n) dt + £ ^(t) cos kt dt = n jn°+2n ^(t) cos ktdt+1 ^(t) cos kt dt = n Jto°+2n ^(t) cos kt dt- Coeficienții Fourier pot fi calculați integrând pe orice interval de lungime 2n an = 1 Jt°°+2n ^(t)cos ntdt „ oricare ar fi t0 E R bn = n Jto°+2n^(t) sinntdt 2 1 15 Seria Fourier asociată unei functii pare este de forma |a0 +22°= an cos nt Seria Fourier asociata unei functii pare este de forma 22° i bn sin nt 2 1 16 Exercițiu (Seria Fourier a functiei “dinț de fierăstrâu” (vezi Fig 2 1)) Sa se arate ca seria Fourier asociata functiei periodice cu perioada 2n si ^>(t) = t pentru t E [—n,n) Serii Fourier 35 3 7 2 7 1 7 -5 5 -1 7 -2 7 -3 7 Figura 2 1: Funcția “dinți de fierăstrău” este Rezolvare Utilizând integrarea prin părți și relația cosn?r = (—l)n obținem an = ± f\t cos nt dt = fZ- t(sin ntj'dt = ± t sin nt\sin nt dt = 0 bn = 7 JCa smntdt = f\t(cosntydt = cos?? t| ^ + iZF JZk cos ?? t dt = — [tt cos îwt - (-%) cos(-n?r)] = cos«7? = (-1)" = (-l)ra 1 2 1 17 MATHEMATICA: FourierTrigSeries [f [t], t, k] In[l]=FourierTrigSeries [t, t, 4] i—> Out[l]=2 Sin[t] —Sin + | Sin —Sin 2 1 18 Se poate arăta (vezi ) că seria (2 3) este convergentă și suma ei este 2 — sin nt = dacă t G Zx ^(t) dacă t Z-zr 0 Funcția p este continuă exceptând punctele tGZ-zr = {fc-zr | fceZ} Se observă că suma seriei Fourier asociate coincide cu tp doar in punctele in care aceasta este continuă Deoarece lim (-1)”-1 - = 0 contribuțiile termenilor devin din ce in ce mai mici pe măsură ce n crește In Fig 2 2 și Fig 2 3 prezentăm sumele parțiale 5 2 10 2 ^(t) = J2(-l)” 1 - sin?? t și S'io(t) = 52(-l)ra 1 -• 36 Aplicații ale transformărilor Fourier și Laplace 3 7 2 7 1 7 -5 ' ' ' ' ' ' 5 -1 ; -2 ; -3 ; Figura 2 2: Suma parțială ^(t) 3 ; -3 Figura 2 3: Suma parțială S'io(t) 2 1 19 Exercițiu Să se arate că seria Fourier asociată funcției tp : [—7F, 7r] —> R, ^(t) = \t\ este prhpcos(2fc+1)t (2-4) Rezolvare Avem «o = I fCr \t\dt = Zf”tdt = TT iar pentru n 0 an = \ t\ cosntdt = f^t cosntdt = [f sin7/t|g — sinntdt = -j-cos ntl =-?-((—l)ra — 1) nz7F |q nz7F w / / bn = țp \t\ sinntdt = 0 (funcție impară) Funcția obținută prelungind R, p(t) = t2 este tj-2 00 ( 1 yn — + 4 52 ——- cos nt (2 5) 3 “ nz n=l Rezolvare Avem a0 = ±JZt2dt=2-ț iar pentru n 0 an = f\t2 cosntdt = f^t2 cosntdt = [t2 sinn t |q — 2 f^t sin ntdt] = Io t sinntdt = [t cosnt|o - f” cosntdt] = ^(-l)ra bn = J22t2 sin nt dt = 0 (funcție impară) Funcția obținută prelungind tp prin periodicitate este o funcție continuă deoarece y?(7r) = y?(—-tf) Din teorema fundamentală (vezi pag ??-??) rezultă că seria (2 5) este convergentă și suma ei coincide cu tp, adică 38 Aplicații ale transformărilor Fourier și Laplace °° f iy« y?(t) = y +4 12 yy- cos nt n= 1 oricare ar fi tGR 2 1 22 MATHEMATICA: FourierTrigSeries [f [t] , t, k] In =FourierTrigSeries [t“2, t, 2] i—> Out[l] = ^-+4(—Cos[t] + iCos ) Figura 2 5: Prelungirea periodică a funcției t* 2 și o sumă Fourier parțială 2 1 23 In cazul in care coeficienții an și bn sunt, numere reale, relația ( tan(/)n = — dacă an cos ntdbn sin nt = Jaf+b'^ cosțnt+(f)n) unde n = — f dacă an = 0 ne permite să scriem seria Fourier asociată unei funcții sub forma |a-o + 12 cosțnt + (/)n) 2 n=0 unde \/an+^n reprezintă amplitudinea, armonicei de ordinul n iar C cu perioada 2% Utilizând formulele gini | ini ini cos nt = - și sin nt = - 2 g 2i obținem Serii Fourier 39 oo oc Z 2ao + E (an cos nt+bn sin nt) = 2 ao + E (an n=1 n=1 ' = 2ao + E (2 (an n=1 ~int i ^-int , ^Ant ~-int e | A e ~e 2 + Vn 2i - ibn)eint +1 (an + ibn)e int) = co + E (cneint+c—ne int n=1 E cneint n= — c unde c0 = 2a0 = 21n J-n ^(t) dt Cn = 2 (an - ibn) = 2n j'-n ^(t) e mtdt c-n = 2 (an + ibn) = j-n ^(t) eintdt adică cn 1 fn = — ^>(t)e—intdt oricare ar fi n E Z 2n J—n 2 1 25 Șirul de funcții periodice cu perioada 2n e—3it e—2it e—it o eit e2it e3it • ••* e , e , e , \j , e , e , e , este ortonormat în raport cu produsul scalar 1 /■n - = ^J ^(t) ^(t) dt- 2 1 26 Definiție Fiecărei funcții periodice : R —> C cu perioada 2n pentru care coeficienȘtii 1 r n Cn — / ^(t) e—int dt 2n J—n există i se asociază seria oc £ Cneint (2 6) n=—c numita seria Fourier complexa a lui ^ 2 1 27 Exercițiu Seria Fourier complexa a functiei “dinti de fierastrau” este £ (-1)n — eint (2 7) n=0 Rezolvare Prin calcul direct, pentru n = 0, obtinem 40 Aplicații ale transformărilor Fourier și Laplace Cn = H-, t e—Mdt = in J-, t dt = ■■ [t e- int |n 1—, e— int 1 L e-inn 2nin e + n einn + in (e—inn ginn 21n [2n cos nn 2 sinnn\ = n - cos nn = -(—1)n nn —n 2 1 28 MATHEMATICA: FourierSeries[f[t], t, k] In =FourierSeries[t, t, 2] i—— Out =ie "—ie" — 1 ie 21t+1 ie21t 2 1 29 Rezultatele prezentate pot fi usor extinse la funcții periodice cu perioada T Șirul de funcții -U, cos wot, sin wot, 2 unde cos2wot, sin2wot, cos3wot, sin3wot, (2 8) este un sistem ortonormat Șsi 2n wo = t în raport cu produsul scalar ■T/2 ^(t) ^(t) dt T/2 ^(x) = 1 k | an = T f-T/2 ^(t) cos nwot dt -ao (an cos nwot+bn sin nwot) C cu perioada T și astfel încât l 0 pentru —T/2 ±oo T nWq contribuțiile termenilor devin din ce in ce mai mici pe măsură ce |?? | crește In Fig 2 8 și Fig 2 9 prezentăm sumele parțiale C i+\ V® 2 sin(»îWo«/2) inwot c i+\ v^lO 2 sin(»îWo«/2) inwot — 2^n=-5 T nuo e 51 — 2^n=-10 T nu0 e in cazul T = 2-tf și a = -tf 2 1 34 Exercițiu (Seria Fourier a funcției triunghiulare periodice (vezi Fig 2 10 )) Fie 0 C cu perioada T și astfel încât 0 1 + - 1 a 3 - f- a o pentru pentru pentru pentru ^(f) = ±oo nliVoaT contribuțiile termenilor devin din ce in ce mai mici pe măsură ce |?? | crește In Fig 2 12 și Fig 2 13 prezentăm sumele parțiale Serii Fourier 45 0 25 0 20 0 15 0 10 0 05 -5 5 Figura 2 11: Coeficienții cn și graficul funcției f în cazul T = a = tt/2 Ss(t) = EL-3 4Sif2Î7aT/2) eira"ot $5(t) = EL-5 4C2y2) eira"ot în cazul T = 2-tf și a = 7r/2 0 8 0 6 0 4 0 2 -5 Figura 2 12: Suma parțială S^t) 2 1 37 Plecând de la orice funcție continuă pe porțiuni f : K b\ —> R putem considera restricția ei la un subinterval [a, /3] C [a, 6], iar apoi putem extinde acestă restricție la R prin periodicitate cu perioada T = (3 — a Seria Fourier corespunzătoare funcției periodice astfel obținute poate fi 46 Aplicații ale transformărilor Fourier și Laplace 0 8 0 6 0 4 0 2 -5 5 Figura 2 13: Suma parțială Ssțt) determinată utilizând formulele prezentate la pag 40-29 In Fig 2 14 prezentăm prelungirea prin periodicitate (perioada T = 3) a restricției funcției f(t) = f2 la intervalul [—1,2] și o sumă parțială a seriei Fourier 2 1 38 MATHEMATICA: In[l] = f [t ] := t*2 alpha = -1; beta = 2; number = 10; T = beta - alpha Plot [f [t] , {t, alpha , beta}, AspectRatio -> 0 7] Show [{Plot [f [t] , {t, alpha , beta}], Plot [f [t - T] , {t, alpha + T, beta + T}] , Plot [f [t + T] , {t, alpha - T, beta - T}]}, PlotRange -> AH , AspectRatio -> 0 3] a[n ] := (2/T) Integrate [f [t] Cos , {t, alpha, beta}] b [n ] := (2/T) Integrate [f [t] Sin , {t, alpha, beta}] Phi [t ] = a /2 + Sum[ a[n] Cos + b [n] Sin , {n, 1, number}] Plot[Phi [t], {t, alpha - T, beta + T}, AspectRatio -> 0 3] 2 1 39 Seria Fourier asociată unei funcții de două variabile ip : R2 —> C periodică cu perioadele (2%, 0) și (0, 2%) ipțajfî) = 3 sau echivalent 1 f2tv f2tv Cn'm = 4^J0 da Jo d(3 C periodică cu perioadele (Ti,0) și (0,72) tp(a, (3) = tp(a+Ti, (3~) = tp(a, (i+To) oricare ar fi (a, /5) G R2 este 48 Aplicații ale transformărilor Fourier și Laplace 2S1 na 2S1 mă e Ti e t2 unde n -^ m=- C cu perioadele T1a1 și T2a2 p(x) = p(x + T1a1) = p(x + T2aa) oricare ar fi x G R2 îi putem asocia funcția periodica p : R2 —> C, țj(x1, x2) = p(x1a1 + x2a2) cu perioadele (T1,0) și (0, T2) Seria Fourier asociata funcției p este Q ■ o • E —2n1 ki ti 2n1 ko to / \ E cU,fc2 e T1 k1T1 e t2 k2T2 (2 12) ki = — *+Tiai) si ei(fe’*) = ei(fc>*+T2a2> este necesar si suficient ca eiqiT1 = 1 si eiq2T2 = 1 adicăa săa aibaă loc relatiile 91T1 G 2nZ si 92T2 G 2nZ care se mai pot scrie 91 G T Z si 92 G TnZ Functia este periodica cu perioadele T1a1 si T2a2 daca si numai daca , 2n , 2n , k G Zb1 + 7TT Zb2 T1 T2 adicăa dacăa si numai dacăa este de forma ^(x1 a1 + x2a2) = e Ti T2 cu k1,k2 G Z 2 1 45 Utilizând notatiile 2n 2n x = X1a1 +X2a2 = (a,£) si R* = Zb1 + — Zb2 T1 T2 seria Fourier (2 12) se poate scrie sub forma kgR* iar formula (2 13) pentru calculul coeficientilor sub forma 50 Aplicații ale transformărilor Fourier și Laplace ■=^1cJL ' unde C este paralelogramul determinat de vectorii T1a1 și T2a2 2 2 Aplicații 2 2 1 2 3 Definitii si rezultate fundamentale 2 3 1 Teoremă Dacă seria 1 2 oo ao + [an cos nx + bn sin nx] n=1 este uniform convergentă și dacă suma ei este funcția f : R —> R atunci oricare ar fi n G {0,1,2, } an = 1 f (x) cos nx dx bn = n J-n f (x) sin nxdx si are loc egalitatea lui Parseval oricare ar fi n G {1,2,3, } (2 15) n 1 °° 1 r 2 a2-y 2+bn) = - yf 2(x) dx- (2 16) n=1 —n Demonstrație Fie sk(x) = ja0 +5Zn=1 [an cosnx + bn sinnx] Deoarece |sk(x) cosnx — f (x) cosnx| = | cosnx| • |sk(x) — f (x)| n —n Serii Fourier 51 Funcția periodică f cu perioada 2n fiind limita unui șir uniform convergent de funcții continue este continua și deci marginită La fel ca mai sus, sk f —^ sk f -2 f2 și R R n n Jf 2(x) dx = lim^c J Sk (x) f (x) dx —n —n — limk^c n n 1 ao I f (x) dx + (an cos nx + bn sin nx)f (x) dx —n —n u™ I" 1,,2i ^k (^2 i l2\1 — 1" 1^21 („2 \ a2\"| — lim^c n 2 ao +2 ^n=i (a„ + b„) — n 2ao + n=i (a„ + b„) zv t 2 ' * n— x\ n nz 2 J n— x \ n n/ 2 3 2 Propoziție Daca f : [—n,n] —> R este o funcție continua pe porțiuni și n tn : [—n,n] —> R, tn(x) — — + Y~^(ak cos kx + Pk sinkx) 2 k—1 un polinom trigonometric de gradul n atunci cea mai mică valoare a integralei hă — [f (x) — t„(x)]2dx (2 17) J —n se obtine în cazul în care ak și Pk sunt coeficienți Fourier (2 15) asociaț funcței f Demonstrațe Utilizand relatiile (??) si (2 15) obtinem A2 — fn f 2(x) dx — 2 fn f (x) t (x) dx + fn t2 (x) dx °n — J—~ f (x) UX 2 J —n f (x) tn(x^x + J — - tn(x^x — f 2(x) dx — aii fn f (x) dx — 2 \ 'n [a/ fn f (x) cos kxdx J n f (x) x a0 j n j (x) o»x Ăd / ^k— 1 ak j n f (x) cos a xUjx + Ș fn f 1%) cin [0/2 + V? Av2 i «2)! I Pk J—n f (x) sin kxdx j + n 2 a0 +2 ^k—1 (ak + Pk ) j — f n f 2(x) dx + n 12(a2 — 2aoao) + Xn—1 (ak + Pk — 2akak — 2Pkbk)] fn f2(x) dx — n [1 a2 i (a2 i b2)1 — j—n J (x) dx n [2a0 k—1(ak + bk +% {2 (ao— ao)2 + z2=1[(ak— ak)2 + (pk— bk)2]} • (2 18) 2 3 3 Teoremă Daca f : [—n,n] —> R este o functie continuă pe portiuni si dacă an, bn sunt coeficienți Fourier asociaț funcței f atunci seria + b^) este convergentă si are loc inegalitatea lui Bessel 1 2 oc a02 i (a2n i b2n n—1 •7T f2(x) dx —n Demonstratie In cazul în care ak — ak si Pk — bk, din (2 17) si (2 18) rezultă 52 Aplicații ale transformărilor Fourier și Laplace Si 1 n 2 ao+y?(afe+ 2 k=i rn ) = [f (x) J —n — tn(x)]2dx > 0 1 oo -a0 ' 12("k + bk) = nlim 2 k=1 1 1 fn 2°0 + + bk) R este o funcție continuă pe porțiuni atunci coeficienții Fourier an, bn asociați funcției f au proprietatea lim an = 0, lim bn = 0 n^c n^c Demonstrație Relațiile rezultă din convergența seriei $2C=1(an + bn) Și din 0 R o funclje continua pe portiuni Dacă modificam valorile pe care le ia f într-un numar finit de puncte din [—n,n] seria Fourier asociata nu se modificaă Făaraă a restrange generalitatea, vom considera doar funcțtii cu proprietatea f (—n) = f (n) Orice astfel de funclje este restrictia la [—n, n] a unei funclji f:R —>R periodice cu perioada 2n (pentru care am pastrat aceeasi notație) 2 3 8 Teoremă Daca f : [—n,n] R este o funcție continuă, derivabilă în intervalele de continuitate și cu derivata f' continuă pe porțiuni atunci seria Fourier asociată lui f este convergenta în orice punct si Serii Fourier 53 1 2 oo ao + [an cos nx + bn sin nx] n=1 f (x—) + f (x+) oricare ar fi x G R 2 Demonstrație A se vedea [?], pag 118 2 3 9 Daca funcția f este continua în punctul x atunci f(x )+/(x+) = f (x) 2 3 10 Teoremă (A doua teorema de aproximare a lui Weierstrass) Orice funcție continua f : R —> R, periodica de perioada 2n este limita unui sir uniform convergent de polinoame trigonometrice Demonstrație A se vedea [?], vol 2, pag 119 2 3 11 Aplicatia r i r -i ?\ a + b b — a V :[—n,n] —> [a,b], V(t) ——^ + -^-t este bijectiva si inversa ei este V-1 : [a,b] —> [—n,n], v—1 (x) = ț—— (2x — a — b) b—a Fiecare funcție f : [a, b] —> R cu f (a) = f (b) se poate prelungi prin periodicitate cu perioada (b — a) pâna la o functie f : R —> R si f = (f ◦ v) ◦ V-1, unde f o V :[—n,n] —> R, (f o V)(t) = + t) este o functie periodica cu perioada 2n Seria Fourier corespunzatoare lui f o v este unde 1 oc -a0 + [an cos nt + bn sin nt] 2 n=1 (2 19) n / \ b a — 1 f f i a+b + b—at| cos nt dx — 2 f f (x) cos nn (2x a b) dx an — n J f \ 2 + 2n U cos — b—a J f (x) cos b—a (Zx a u)Ujx —ny/ a n / x b b — 1 f f i a+b i b—at| sin nt dx — 2 f f (x) sin nn (2x a b)dx Un — n J f \ 2 + 2n U sin — b—a J f (x) sin b—a (^x a u)Ujx —ny/ a Deoarece f — (f ov) ov—1 efectuandân (2 19) substituia t — bna(2x — a — b) obtinem seria Fourier corespunzatoare lui f nn , , nn cos (2x — a — b) + bn sin - b — a b — a (t) o cn în loc de coeficienții Fourier ai funcției ^>(t) sunt cn 2 3 13 Teoremă Avem ^(t) o (pn fi(t) O fin Demonstrație Afirmația rezultă din relația /•T/2 /-T/2 / (a^(t) + fifi(t))e inw°tdt = a / ^>(t) e J-T/2 J-T/2 aip(t) + fifi(t) o a^n + fifin r T/2 inw°t dt + fi fi(t)e-inu°tdt J-T/2 2 3 14 Teoremă Avem ^(t) O C-n F(t) O Cn ^(t-tfi) o e inW0t°cn F(~t) O c-n 54 Capitolul 3 Transformarea Fourier discreta 3 1 Definiție si exemple 3 1 1 Definiție Spunem despre o funcție de variabilă discretă : Z —H C că este periodica cu perioada N dacă ^>(n+N) = ^>(n) oricare ar fi n E Z 3 1 2 Oricare ar fi k E{0,1, ,N—1}, funcția exponențiala Z —H C : n H e ~ kn este o functie periodica cu perioada N deoarece 2ni k(n I IN) i kn 2kni 2ni kn e n k(niN) = e n kn e2 = e n kn Am utilizat relația e2kni = cos 2kn + i sin 2kn = 1 3 1 3 Orice functie : {0,1, ,N—1} —H C se poate prelungi prin periodicitate pana la o functie periodica : Z —H C 55 56 Aplicații ale transformărilor Fourier și Laplace cu perioada N și orice funcție periodică : Z —> C cu perioada N este complet determinata de restricția ei la mulțimea {0,1, ,N — 1}, adica de numerele ^(0), ^(1), ^(N—1) Deoarece = (n) = C este funcție periodică cu perioada N atunci , , 1 2ni „k N J 2ni , d(n) = N e n nk e n kmip(m) k=0 m=0 / x 1 1 2ri „fcN—1 2ri z x C cu perioada T f (t+T) = f (t) oricare ar fi t G R neteda pe portiuni este dezvoltabilă în serie Fourier oo 1 z t n kt 1 a 2ni kt + (Q >1 f (t) = cke T cu ck -p e T f (t) ut (3 2) T0 k= o Din (3 1) rezultă ca :Z —>C periodica cu perioada N admite reprezentarea Transformarea Fourier discreta 57 TV-l N-1 / \ X—2ni nk 1 X—2ni km / \ T(n) — 52 cke N cu ck = n 52 e N T(m)- (3 3) k=0 m=0 3 1 6 Definiție Transformata Fourier a funcției periodice T : Z —C cu perioada N este funcția TV-l F[t] : Z —> C, F[ C este functia F[t] : {0,1} -^ C, F[T](k) — e-knT(n) — (-1)knT(n) n=0 n=0 adica avem F [T](0) — T(0) + T(1) F [T](1) — T(0) - T(1) 58 Aplicații ale transformărilor Fourier și Laplace 3 1 10 In cazul N = 3 avem 1 daca kn G 3Z — 1 daca kn G 3Z+1 — i + i^3 daca kn G 3Z+1 și prin urmare transformata Fourier a unei funcții : {0,1, 2} —> C este functia F[ ] : {0,1, 2} —> C, 2 F[^](k) = e-kra^(n) n=0 adica avem F M(0) = ^(0) + ^(1) + ^(2) FM(1) = ^(0) + (—1 — P33) ^(1) + (—1 + p2t) ^(2) F m(2) = ^(0) + (—1+i^) ^(1) — (—1+^(2)- Ultimele relații se mai pot scrie / F[^](0) \ F[^](1) I \ F[^](2) J 1 1 1 \ / ^(0) \ 1 — 1 — — 1+^23 l I • , 1 — 1 +i — 1 — i I \ w(2) / \ ± 2 + 1 2 2 1 2/ ^vy 3 1 11 In cazul N = 4, transformata Fourier a unei functii : {0,1, 2, 3} —> C este functia 3 3 FM : {0,1, 2, 3} > C, FM(k) = y>-kn^(n) ^^(—i)fcra^(n) n=0 n=0 adica avem F M(0) = ^(0) + ^(1) + ^(2) + ^(3) F M(1) = ^(0) — M1) — ^(2) + M3) F M(2) = ^(0) — ^(1) + ^(2) — ^(3) F M(3) = ^(0) + M1) — ^(2) — M3) Transformarea Fourier discreta 59 3 1 12 In relațiile (3 4) și (3 5) putem considera că n parcurge mulțimea Zn a claselor de resturi modulo N Astfel, transformata Fourier a funcției T : Zn —C este functia F[t] : ZN C, FM(fc)= £ e-2^ kn^(n) nEZN iar transformata Fourier inversa a lui t este functia F-1[t] : Zn -^ C, F-1[^](fc) = N £ e& kn^(n) nEZN De exemplu, în cazul în care N = 2M + 1 este impar, multimea {-M, -M + 1, ,M- 1,M} este sistem de reprezentanti pentru Zn si prin urmare F[ C cu perioada N se poate identifica cu CN asociind funcției p elementul (xq, xn-i) = ( CN : x = țx0,xi , ,xN 1') F[x] = țF[x]0, F[x]!,F^n-^ unde w-i F\x]k = n=0 Transformarea inversă este w-i F~1lx]k = -^^knxn n=0 Transformarea Fourier discretă 61 3 1 15 Exercițiu Fie mG{0,1, , N — 1} fixat și funcția delta discretă ăm : {0,1,N — 1} —> C, 1 dacă ii ni 0 dacă ii ■ in Să se determine F Rezolvare Avem w-i F\5n0]W = E e-^1 kn6m(n) = e~^km n=0 In particular, scriind 6 in loc de 5o, avem F (fc) = 1, F (fc)=e-^fc, F (fc) =e~^k, etc Figura 3 2: Funcția 63, partea reală (stânga) și partea imaginară a transformatei Fourier Ffâs] in cazul 7V = 10 3 1 16 Exercițiu Să se calculeze transformata Fourier a funcției tp : {0,1, , N — 1} —> C, C, ^(n) = CN—i- Sa se calculeze F[^] Rezolvare Utilizand formula (binomul lui Newton) m (a + b)m = ^2 Cm am—nbn n=0 obtinem N - 1 y F[ R este parăa Un) = ^(n) atunci, utilizând schimbarea n —n, obtinem F [^](fc) = E e— 2N kn^(n)= £ e 2N1 kn^(—n) n=—M n=—M = £ e2-N kn^(n) = F[^](fc) n=—M adică transformata Fourier F[ C este funcția F[t] : {0,1} —> C, F [T](0) = T(0) + T(1) F[t](1) = t(0) - t(1) 3 2 2 In cazul în care perioada este un numar par 2N, unei functii T : {0,1, , 2N-1} —> C îi putem asocia functiile Tpar : {0, 1Î •••, N — 1} -> C, Tpar (n) = T(2n) T impar : {0, 1j •••, N — 1} -C- Timpar (n) = T(2n + 1) si avem 2N — 1 2ni i F[T](k) = E e-ÎN knT(n) n=0 — e- 2N k 2mT(2m) + V'' e— 2N k (2m+1) t(2m + 1) — / j e Ti^mj ~+ / j e Ty^’m i- m=0 m=0 N—1 — — km / x i — ni k N—1 — — km ( x = Y e N km Tpar (m)+e N ^2 e N kmTimpar (m) m=0 m=0 J F[Tpar](k) + e — Nk F[Timpar](k) dacă k G{0,1, ,N—1} F[Tpar](k — N)+e N k F[Timpar](k — N) daca kG {N, N+1, , 2N — 1} Utilizand formula obȘinuta, calculul transformatelor Fourier ale funcțiilor de forma T : {0,1, , 2n — 1} —> C se poate reduce succesiv pana la calculul transformatelor unor functii de forma : {0,1} —> C Aceasta metoda de calcul foarte eficienta este numita transformarea Fourier rapida 64 Aplicații ale transformărilor Fourier și Laplace 3 2 3 Calculul transformatei Fourier a unei funcții : {0,1, 2, 3} —> C se poate reduce la calculul transformatelor Fourier ale functiilor Fpar : {0, 1} -> C- ^impar : {0, 1} -> C Fpar (n) = ^(2n) ^impar (n) = ^(2n + 1) Avem f F[^par](k)+e 2 k F[^Pimpar](k) daca k G{0, 1} FM(k) = { ni k [ F[Ppar](k-2)+e 2 kF[Pimpar](k-2) daca kG{2,3} ^>(0) + ^>(1) + ^>(2) + ^>(3) dacă k = 0 ^>(0) — i^(1) — ^(2) + i^(3) daca k = 1 ^>(0) — ^>(1) + ^>(2) — ^>(3) dacă k = 2 ^>(0) + i^>(1) — ^>(2) — i^(3) daca k = 3 3 2 4 In cazul unei funcții de forma : {0,1, 2, , 23 — 1} —> C functiilor Fpar : {0,1, 2, 3} —> C, ippar (n) = >p(2n) ^impar : {0, 1 2, 3} C- ^impar (n) = lf(2n + 1) li se poate aplica direct formula obfinuta la punctul anterior 3 3 Proprietăți ale transformării Fourier 3 3 1 3 4 Funcții proprii ale transformării Fourier 3 4 1 Pe parcursul acestei secfiuni vom utiliza pentru transformarea Fourier definita F[^l(k) — —1— e n kn^(n) cu F 1[^](k) — —1— 5^ e n kn^(n) 1 ) Z-i e N w(n) cu 1 ) Z-i e N w(n)- VN n=0 VN n=0 Transformarea Fourier discreta 65 3 5 Transformarea Fourier bidimensională 3 5 1 Fie N1,N2 E{2, 3, 4, } Dacă k1, k2 EZ atunci funcția 2ni kn-i 2n1 kirUn : Z x Z —y C, ^(ni,«2)=eNi eN2 2 2 (3 8) este o functie periodica cu perioadele (N1,0) și (0, N2), adica ^(«i,«2)= > k uJ -D'i —i N2 —i 2ni k im 2ni C, F[T](k1 ,k2) = Ț? e-Nîkînîe-N2 k2n2T(n1,n2) nî=0 n2=0 Transformata Fourier inversa, a lui t este funcția I Nî-1N2 — 1 2 2 F—1[t] : ZxZ C, F-1[^](k1 ,k2) = —— £ e N 1 îeN2 2 2^,«2) N1N2 nî=0 n2=0 3 5 5 Transformatele Fourier ale funcțiilor periodice cu perioadele (N1,0) si (0, N2) sunt functii periodice cu aceleasi perioade si F [f—1[t]] = t si F-1[F [t]] = t- Astfel, transformata Fourier a funcției T : Z Nî x Z N2 —C este functia F[t] : ZNî x ZN2 —> C definita prin relatia , ,1 i/, , \ -Ip kînî - fe2«2 / \ F[T](k1,k2)= > y X e Nî e N2 T(n1,n2) nî^Z^î n2€ZN2 iar transformata Fourier inversa functia F—1[t] : Znî x Zn2 —> C, F [T](k1 ,k2 ) = NN X X e Nî e N2 T(n1,n2) 1 2 nî€ZNî n2 €Zn2 3 5 6 Relatiile din teorema prezentata la pag 65-3 permit reobtinerea exacta sau aproximativă a unei functii din transformata ei Fourier In figura ?? prezentam în diverse reprezentari grafice functiile : {-3, -2,-1, 0,1, 2, 3} x {-4, -3, -2,-1, 0,1, 2, 3, 4} —> R definite prin f (n1,«2) |nm2| «2 + «2 +1 și Transformarea Fourier discretă 67 2 2 ^(7ii,7i2) = î5 E E k±=—2 /c2 =—2 7711= —3 7712 = —4 3 5 7 MATHEMATICA: 68 Aplicații ale transformărilor Fourier și Laplace In = M1 = 3; M2 = 4; N1 = 2 M1 + 1; N2 = 2 M2 + 1 f[n , k ] := N[Abs[n k]/(n*2 + k~2 + 1)] Ff[n , k ] := N[(1/(N1 N2)) Sum[Exp Exp f[a, b], {a,-M1, M1}, {b,-M2, M2}]] g[n , k ] := N[Re[Sum[Ff[a, b] Exp[-2 Pi I a n/N1] Exp[-2 Pi I b k/N2], {a,-2,2}, {b,-2, 2}]]] DiscretePlot3D[f[n, k], {n, -M1, M1}, {k, -M2, M2}] DiscretePlot3D[g[n, k], {n, -M1, M1}, {k, -M2, M2}] DiscretePlot3D[f[n, k], {n, -M1, M1}, {k, -M2, M2}, ExtentSize -> Full] DiscretePlot3D[g[n, k], {n, -M1, M1}, {k, -M2, M2}, ExtentSize -> Full] Image[Table[f[n, k], {n, -M1, M1}, {k, -M2, M2}] ] Image[Table[g[n, k], {n, -M1, M1}, {k, -M2, M2}] ] 3 5 8 Fie a1,a2 GR2 doi vectori liniar independenți, Ni,N2 g{2, 3, } și R = Za1 +Za2 = { n1a1 + n2a2 | n1,n2 G Z } rețeaua Bravais cu baza {a1,a2} O funcție periodică : R —y C cu perioadele N1a1 si N2a2 poate fi considerată ca fiind o funcție : R —y C definita pe mulțimea cu N1N2 elemente R = ZN1 a1 + Zn2 a2- Functiei îi putem asocia functia periodică : ZN1 x ZN2 —y C, C cu perioada T f (t+T) = f (t) oricare ar fi t G R neteda pe portiuni este dezvoltabila în serie Fourier oo f(t)= Y Cn ent (3 11) n=—c cu 1 rT/2 Q ■ 1 cT „ c =A / f(t)e— nt dt = — / f(t)e—nt dt Cn = rri f (t)e 1 dt = rri f (t) e 1 dt T J—T/2 T Jo Transformarea Fourier discreta 71 3 7 2 Relația (3 11) se mai poate scrie 1 OO T 1 V e — nt f(f)e- nt dt = f(t) (312) f y e T / J (t) e T dt — J (t) (3 12) T n=-O 3 7 3 In aplicații, funcția f este determinata “experimental” si se cunosc valorile doar in anumite puncte t G Admitand ca stim doar valorile f(0), f (T) , f (2N , - f ((N-1)N unde N este un întreg pozitiv fixat, putem aproxima integrala T f (t) e-2r nt dt 0 cu suma Riemann N-1 /, x T E ((k+1)N k=0 v k w - nkT T E-1 - — nk e t nk n — N e n nk k=0 Avem Cn — limN -o N 5- e " nk k=0 dar pentru N fixat, aproximarea 1 N-l z rp 1 2Z1 nk e ( 7 T Cn — > e N nk f k — nN N k=0 ar conduce la 1 N 1 / C«+N « £ £ e- 2N k f (k k=0 în contradicție cu N— 1 z rp \ nkf YN) ~ Cn k=0 lim cn — 0 n—±o 3 7 4 Utilizand formula (suma seriei geometrice) N N-1 Y qm — 1 + q + q2 + ■■■ + qN-1 m=0 daca q — 1 1-qN 1-q dacă q — 1 obtinem însa relatia -i N—1 N—1 z rp \ z rp 1 2ni nm mk e ( j T \ e d v > e2f (d-1)- , este o matrice unitarăa deoarece d-1 d-1 1 -ni kn — —nL mn 1 221 (k — m)n c e d kne d mn = d}^ e d (k m)n = $km = d n=o d n=o dacaă dacăa k = m (modulo d) k = m (modulo d) 1 0 Transformarea Fourier discreta 73 3 7 7 In cazul unui spațiu Hilbert d-dimensional este convenabil uneori sa indexăm coordonatele folosind inelul Zd al claselor de resturi modulo d Deoarece 2ni kn 2ZÎ (k+id)n k(n + îd) e d kn = e d (k+jd)n = e d k(n+jd) oricare ar fi j E Z, au sens relațiile 1 V—2ni kn 1 X—— 2ni kn yk=^E ed xn, xn=^E e d yk■ V d neid Vd keZd 3 7 8 Orice spațiu Hilbert d-dimensional V poate fi identificat cu Cd alegand o bază ortonormata {|vn\}neZd Spatiul V poate fi astfel identificat cu spațiul funcțiilor de forma : Zd —— C Transformarile Fourier directa si inversa se pot scrie F = — E e^MK|, F-1 = — £ e-Tte|t^| = F*■ V d k,nezd V d k,nezd 3 7 9 Deoarece F4 = I valorile proprii ale lui F apartin multimii {1, -1, i, —i} 3 7 10* Se stie (a se vedea pag 86-13) ca pentru transformarea Fourier continua Zoc e1?x^(x) dx functia x2 fk : R —— R, fk(x) = Hk(x) e 2 definita cu ajutorul polinomului Hermite Hk, este functie proprie 1 Z 2 £2 dx e1?x Hk (x)e- V = ik Hk «)e-~ 3 7 11* Folosind functia periodicăa Fk : R —— R, Fk (x) = E hE (ad+x)) Q=-^ ' V ' z t x 2 -1 ( Vd (ad+x) ) e2 d cu perioada d definim functia Zd —y R : n i—y / t \ 2 Fk(n)= g HJ(ad+«)} e-H{ad+n) a=- c r(a+1)V2nd e a— c r (a+l^/ 2nd a— R cel mult d pot fi liniar independente 3 7 14* Funcția lui Jacobi 03 definită prin relația [?, ?] oc 03(z,T) = 52 einTa2 e2niaz, Im (t) > 0 a=—c este importanta pentru fizica teoretica datorita proprietatilor ei, dintre care menționam 03(z + m + nT,T) = e—inTn2 e—2ninz 03(z, t) o £ ■ \ 1 ——A (z i A 03(z, iT) = —= exp T 03 — VT \it r 3 7 15* Oricare ar fi k G (0, to), se poate arata [?, ?] că 03 k i«ă dd Kd d—1 kn 03 n=0 1 n d Kd Din aceasta relație rezulta ca functiile (vezi Figura 3 4) oc gK (n) e— (ad+n)2 a=-c ni (F Kd gK : Zd > Rî Figura 3 4: The functions gi/3, gi și g3 in cazul d = 31 Figura 3 5: Funcțiile e s-1 , e 2 care pot fi privite ca versiuni discrete ale funcțiilor (vezi Figura 3 5 ) pK:R—>R, pK(x)=e 2^ kG(0, 00) verifică relația ^[g«] = -Ugi y K K 3 7 16 Dacă alegem pentru transformata Fourier a unei funcții continue definiția atunci (a se vedea pag 80-9) 76 Capitolul 4 Transformarea Fourier a functiilor 4 1 Exemple 4 1 1 Definiție Fie : R —> C Funcția Zoc (în cazul în care există) se numește transformata Fourier a lui None] Plot[FourierTransform[HeavisidePi[t], t, x], {x, -30, 30}, PlotRange -> All] 4 1 3 Exercițiu Fie a G (0, to) și : R > R tot) = J 1 daca |x| a * : R R, ^(x) = j o daca |x| > a (4 1) Sa se arate că (vezi Fig 4 1) Ft^]( Out[l] = J=Sinc [^1 V 2tt L z j Figura 4 2: Funcția (4 2) in cazul a = l și transformata ei Fourier 4 1 5 Exercițiu Fie a G (0, oo) și ^(x) = 1 -0 Să se arate că (vezi Fig 4 2) = LALLL dacă |x| a (4-2) eit e ii COS t = Rezolvare In cazul 0, utilizând schimbarea de variabilă x f 9 t 1 — cost sin2 - = 2 2 —x și formulele obținem ^M«) = JX e^țx) dx = J°a ete« (1 + dx + foa eia* (1 - dx = fo (1 - □ d x + Joa eij'e (1 - □ dx = 2 /“ (1 - □ cos xț dx = (sinx^'dx= | (1 - ^) sin^|“ + J J“ sin^drr = cosași" = ^(1 -cos^) = 4sing/2) In cazul £ = 0 avem ^M(O) = J°a (1 + 1) dx + foa (1 - l) dx = a Transformarea Fourier a funcțiilor 79 Deoarece lim 4 sin2 («4/2) «42 sin(a4/2)l2 = a lim L «4/2 = a transformata Fourier este o funcție continuă 4 1 6 MATHEMATICA: Definiția utilizată eitx V 2tf J —OO In :=FourierTransform[(l-Abs[t]) HeavisidePi[t/2], t, x] i—> Out[l]= ~3~I~Ș ,'— V2tt xz Figura 4 3: Funcția e 1*1 și transformata ei Fourier 4 1 7 Exercițiu Să se arate că 2a a2 + ț2 oricare ar fi a G (0, oo) Rezolva,re Considerând integrala in sensul valorii principale avem jr[e-«l-d] (£) = dx = e-«l-d (cos țx | j sjn ^x = fZ e “l:l’l cos £x dx = 2 Jo°° e~ax cos £x dx Integrând de două ori prin părți obținem relația e~ax cos £x dx = |e “* sin4a?| +| Jo°° e “‘T sin țx dx IOO 2 QQ , 2 00 = — ^e~ax cos Jo e~ax cos dx = ■£? — țț Jq e~ax cos Sțx dx adică [ e~ax cos £x dx = — y-r I e~ax cos £x dx , Jo ț2 C2 Jo din care deducem o e ax cos ^x dx = a2 + Out[l]=e AbsH In :=FourierTransform[Exp[-Abs[t]], t, x] Figura 4 4: Funcțiile e T"/2, e e 2,'~ (stânga) și transformatele lor (dreapta) 4 1 9 Exercițiu Să se arate că oricare ar fi a G (0, oo) Rezolvare Avem e-« r2ei 2 Plecând de la integrala e 2 e 2 e~az2dz = 0 a funcției fțz) = e “~2 de-a lungul drumului dreptunghiular din Figura 4 5 arai am câ Transformarea Fourier a funcțiilor 81 Avem /*r /* oo li^ / e—at dt = e~ai? dt r^oJ-r J—o Alegând pentru drumul liniar ce unește r cu r — i2^ parametrizarea £ Y1 : —> C, Yi(t) = r — it obținem relatia rr— „ c 1 z t x2 £ £ Q C1 t2t2 / e—2dz = e—a(r—it~) (—i)^dt = — i^ e—ar eirt?+~£rdt Jr Jo 2a 2a Jq din care rezulta lim / 2a e—az2 dz = 0 r^°J r —r r ± i 2a1 Figura 4 5: Drumul dreptunghiular utilizat Similar se arata ca —r 2 az2dz = 0 —r— 2a Alegand pentru drumul liniar ce uneste —r i 2-a cu r i 2« parametrizarea Y2 : [—r, r] —> C, Y2(t) = t i£ 2a r^o e obținem relatia r r— 2a az2 dz = dt i 2 a e —r— r din care rezulta o r^o r— 2a az2 dt —r— 2a —o e e ra{t—i 2a) 82 Aplicații ale transformărilor Fourier și Laplace 4 1 10 MATHEMATICA: Definiția utilizata F [^](x)=■—= eitx (u) du] df = — ^(x) -tt \«/ — tt / k și ftt / ftt I eiKW / e— — tt \J—tt iKUX^(x) dx) du = —p(f) ) k Demonstrație Oricare ar fi a G (0, w) avem \ e—“«2—iK?x (JX eiK«>(u) du) df = ftt m(u) 17tt eiKUu—x) e—“£2 df] du = 7^ ftt m(u)e— K ' du = J—^ ^(u) [J — tt e e “Sj uu = V a J — tt ^(u)e ' uu' Utilizând în ultima integrala schimbarea de variabila u = x + 2^t obținem relația Ztt 9 / r \ - /n r ( ra \ 2 e—a^ —iK?^ eiK?u^(u) du] df = U x + -^t] e—t dt -tt \J—tt J k J—^ y k J care pentru a \ 0 devine D e—^( D eiK?u^(u) df df = D ^(x)e—t9 dt J—tt \7—tt / k J—tt Dar V? ftt t2 2Jn ftt P 1 2n / ^>(x) e t dt = ^>(x) / e t dt = — ^>(x) k J—tt k J—tt k A doua relatie din enunțul teoremei se poate demonstra similar Transformarea Fourier a funcțiilor 83 4 3 2 Definițiile transformărilor Fourier directă și inversă se bazeză pe teorema precedentă Există mai multe posibilități de a defini aceste transformări, trecerea de la o variantă la alta făcăndu-se foarte usor 4 3 3 In cazul k = 1, relatiile din teoremă se pot scrie f°0o e i?x (j 0o ei?>(u)df = ^(x) JX ei?u J °0o ■ dX du = ^(f) si sugerează pentru transformările Fourier directa și inversă alegerile F[^], F[^](f) = J 0o ei?x ^(x) dx F-1 [ț], F 1[ț](x) = f 0o e i?x ^(f) de Ea este alegerea preferată în unele cărti de ecuatiile fizicii matematice si va fi utilizată pe parcursul acestui capitol, exceptând ultima sectiune dedicată aplicaților în mecanica cuantică 4 3 4 In cazul k = 1/ă, relatiile din teoremă scrise sub forma — 1 e i^x/^ [ —1 fo ei^u/Rm(u) d''u\ df — m(x) x2:: ' o e U 2n ■' o e ^(d) ' J df = ^(X) —1— i° ei^u/^ f —1—fi 1° e iux/Rw(r) dx^ du — m(f) ț2nă I o e \F2nh 2n J o e (u) d^ df = (x) dx aparține de asemenea spațiului S(R) si au loc relațiile (F M)(fc) = F [(ix)k F [^(fc)] = (-i C este o funcție indefinit deriv-abila si cu derivatele functii marginite Relatia F[T(k)](f) = / (x) dx = (-if)M x dx = (-if)kFM(£) J—ro J— obtinuta utilizand integrarea prin părți conduce la egalitatea |fkf[V](î) | = |f'V i(«)| care arată ca F[t] E S (R) 4 3 9 Teoremă Transformarea Fourier a funcțiilor de proba Zoc e'ej"p(x) dx este o aplicație bijectiva si inversa ei este transformarea 1 c F—1 : S(R) -^ S(R) : F—1 [^], F—1 [^](x) = — df 2n 7 — ^ Demonstrare Afirmata este o consecința directa a teoremelor anterioare 4 3 10 Se poate arata ca transformările F±1:S(R) —>S(R) sunt continue, adica Tn T =^ F±1[^n] ->F±1[^] 4 3 11 Transformarile Fourier directă si inversă au expresii foarte asemanatoare Utilizand schimbarea de variabila f = — y obtinem F—1[^](x) = - e—^(f) df = - eixy^(-y) dy = -F[ C, fi(y) = fi(-y) In particular, F [F [v]](p) = 4 3 12 Polinoamele Hermite (a se vedea pag ??-??) 2 dn 2 Hn(x) = (-1)n ■ e , n G {0,1, 2, } verifica relațiile de recurența Hra+i(x) - 2xHn(x) + 2nHra i(x) = 0, Hn(x) = 2nHra i(x) 4 3 13* Teoremă Oricare ar fi nG{0, 1,2, }, funcția x 2 in : R —> R, 1n(x) = Hn(x) e 2 este o funcție proprie a transformării Fourier F x2 - ț2 Hn(x)e ~ (fi) = ypîn in Hn(f)e ~ Demonstrație Din relația (F[v])(k) = F[(ix)kv] care se poate scrie dk F [xk v](f) = (-i)k F |v|«) rezultă F Hn(x)e x2 2 = Hn 1 F e x2 2 ți prin urmare (a se vedea pag 84-6) F x2 Hn(x) e- 2 = \ 2n //„ -i^ e 2 \ df/ Folosind metoda inducției matematice vom arata ca / d \ ț2 ț2 Hn (- df) e ' =in Hn«)e 2 Relația are loc pentru n = 0 si presupunand ca / d \ ț2 ț2 Hk ț-idfj e 2 =ikHk(f)e 2 pentru orice k ; d TT I : d\ ,x— •)/ ) I l \ / / l i d \ „ — yy — -2i d?Hn— 1 i d?J e 2 - 2(n- 1) Hn—2 i d? J e 2 £2 - in—1 Hn— i(€)e— - £2 - 2(n- 1)in—2 Hn—2«)e— t j2 = in l-ZH'n—i® + 2(Hn—1(0 + 2(n-1) Hn—2 (0] e— “ 2i d? i’ d? = in [-4(n-1) Hn—2(() +2(Hn—1(() + 2(n-1) Hn—2 (()] e— V £2 £2 = in e— t = in Hn(() e— “ 87 88 Aplicații ale transformărilor Fourier și Laplace Capitolul 5 Transformarea Fourier a distribuțiilor 5 1 Exemple 5 1 1 Se știe că y/—4 nu există în R Dacă scufundăm pe R în C R C : x x + 0i identificănd numărul real x cu numărul complex x + 0i atunci 4 există, dar nu este un număr real Similar, o funcție nederivabilă clasic poate deveni derivabilă dacă scufundăm funcțiile uzuale (doar o parte dintre ele !) în spațiul distribuțiilor, dar derivata nu mai este o functie uzuală ci o distribuție 5 1 2 O cunoastere în profunzime a teoriei distributiilor este ceva mai dificil de realizat, dar pentru a efectua anumite calcule cu distributii este suficient săa cunoațstem cateva elemente simple: • O distribuție (temperată) este o functie liniară si continuă f : S(R) —> C : (se preferă să se scrie în loc de f(^)) definită pe spatiul vectorial S (R) format din funcțiile indefinit derivabile : R —> C : x ^>(x) 89 90 Aplicații ale transformărilor Fourier și Laplace care, atât ele cât și derivatele de orice ordin, descresc la ±w mai rapid decât 1/|x|n, oricare ar fi nEN Distribuții de tip funcție: distribuia corespunzatoare unei functii f : R —> C în cazul în care exista, este Tf : S(R) —> C, , P) d= [ f (x) P(x) dx J —tt (in mod uzual, distribuia corespunzatoare lui f se noteaza tot cu f) Distribuția Dirac: pentru fiecare a E R, aplicatia da : S(R) —> C, {dapp} = p(a) este distribuie In loc de do se scrie doar d, adica avem (d, p} = p(0) Derivata distributiei f : S(R) —> C este distribuia f' : S(R) -^ C, f» d=f — (f,p‘) Produsul dintre funcția indefinit derivabila d : R —> C si distribuia f : S(R) —> C, în cazul în care exista, este distribuia df : S(R) C, (df,p) d=f (f,dp) Transformata Fourier a unei functii din S(R) este tot o functie din S(R) si F : S (R) —> S (R) : p F[p] este transformare bijectiva Transformata Fourier a distributiei f este F[f ] : S(R) -^ C, (F[f ],p) d=f f F[p]> iar transformata Fourier inversa a distributiei f este distribuia F—1[f] : S(R) -^ C, F‘ 1 f d=f f F‘ 1 , 5 1 3 Se poate arata ca functia (vezi Fig 5 1) p : R —> R, p(x) = apartine spatiului S (R) Avem p'(x) = -2x e—x2, p'! (x) = (4x2 x2 e 2)e-x2 Transformarea Fourier a distribuțiilor 91 Figura 5 1: Funcția p(x) = e și in general, p^ (a?) = polinom ■ e Știm că oricare ar fi xgR înlocuind x cu x2 obținem oricare ar fi x G R Oricare ar fi n G {1,2,3, }, avem 5 1 4 Exercițiu Să se arate că relațiile Tf : 5(R) —> C, Sa : ă>(R) —> C, pentru x —> oo = f f(x)p(x)dx (6a,p) = p(a) definesc aplicații liniare Rezolvare Avem {Tf,ap + (3f) = f f(x)(ap(x) + 0^(x)) dx = a(Tf, p) +ț3(Tf,^) (6a, ap + (3ip) = țap + /3 f)ța) = ap(a) + Șipța) = a{6a, p) + (3{6a, tp) 5 1 5 O funcție este numită funcție derivabilă dacă este derivabilă in fiecare punct Funcția Heaviside (vezi Fig 5 2) 92 Aplicații ale transformărilor Fourier și Laplace 1 r H(x) Figura 5 2: Funcția Heaviside H H : R —> R, H (x) = 0 daca x 0 nu este o funcție derivabila deoarece nu este derivabila în 0 h/( ) I 0 daca x = 0 x | nu exista daca x = 0 Distribuția corespunzatoare (numita distribuția Heaviside) H(x) ^>(x) dx = / ^>(x) dx Jo este însa derivabila si (/h= ~(H, 0 Să se arate ca (Tf y = Tf, +25 Rezolvare Integrând prin păr(i obtinem ((Tf y, + = -(Tf, + = - JX f (x) +'(x) dx = î'"x x2 +'(x) dx + f0^(^x + 2)+'(x) dx = -x2 +(x) |-^ + 2x +(x) dx -(^x + 2)+(x)|~ + J(T 277+(x) dx = 2+(0) + JX f'(x) +(x) dx = {Tf + 25, TI unde f' este derivata clasica {2x daca nu exista dacă —daca x 0 prelungita arbitrar în x = 0 5 1 7 Exercițiu Știm că (sin x) = cos x Să se arate ca derivata distributiei corespunzatoare functiei sin x este distribuia corespunztoare functiei cos x, adica 94 Aplicații ale transformărilor Fourier și Laplace (Tsin x ) — Tcos x ■ Rezolvare Integrând prin părți obținem ((Tsin x) , p} — (Tsin x, p } — J—o sin x p (x) dx — — sin xp(x)|°o + JO cos xp(x) dx Deoarece lim p(x) — 0 x^±o obținem Zo cos xp(x) dx — (TCosx,p) o oriacreâ ar fi p G S (R) si prin urmare (Tsin x ) — Tcos x ■ Se observă câ p a jucat doar rolul unui catalizator, dispârând din rezultatul final 5 1 8 Exercițiu Sâ se arate câ dacâ d G Co(R) atunci dda — d(a) da Rezolvare Avem (dda, p) — (da, dp) — (dp)(a) — d(a) p(a) — d(a)(da, P) — = (1, F[t]) = Out|l|\/2nD1racDelta[x] 5 1 14 Exercițiu Sa se arate ca F 2 (Sa + S-a) =cos ax, F [Sa]=e1ax, F [cos ax] = n(Sa + S-a) Rezolvare Avem 96 Aplicații ale transformărilor Fourier și Laplace ri 2 eiaxp(x) dx = (elax ,p) -o Transformarea Fourier fiind liniară obținem F [|(do + d-a)] =1(F[ă„]+ F[M) = Din relația precedenta rezulta ca F -1[cos ax] = |(ăa + d-a) Dar efectuând schimbarea de variabila x — —x obținem (F[cos ax], p} = (cosax,F[p]) = J»cosax (J» elxtp(t)dt^ dx elax | e-lax = cos ax 2 = fcos ax e lxtp(t)dt^ dx = (cos ax, 2nF 1 [p]) = (2nF 1[cos ax],p) adică F [cos ax] = 2nF 1 [cos ax] 5 1 15 MATHEMATICA: Definiția utilizata F [^](x) = ■—= eitx o converge la distribuia f daca lim n^o fn = f lim (fn,p) = (f,p), oricare ar fi pgD(R) Transformarea Fourier a distribuțiilor 97 Figura 5 4: Graficele funcțiilor /i, fa, și fa 5 1 18 Exercițiu Funcției (a se vedea Figura 5 4) fn : R — [ Ș dacă kr ii corespunde distribuția Tfn : P(R) C, /•oo n fl/n fafa ’fa = ffafa fafa dx = - fax) dx J—oo 4 J—l/n oricare ar fi n G N* Limita șirului de funcții (Jn)n>i este f ' M IR, f oc dacă x = 0 = = diu:ă x^0 dar lim Tf =5 Jn Rezolva,re Utilizând schimbarea de variabilă t = nx obținem limn >oo(7/n, fa = lim^oo Ș fax) dx = | lim^oo j]x tp (±) dt = | j]x R, fn(x) = n x îi corespunde distribuția regulată / x > 1 fx sin nx , Tf ■ D(R) C, (Tfn , = - —— ^(x) dx n J — X x oricare ar fi n G N* Limita sirului de funcții (fn)n>1 z „ f w daca x = 0 f : R —h R, f (x) = lim fn(x) = C : dx -x x nu defineste o distribuie Se poate însa arata ca functionala P^: D(R) —H C, (r1 = lim ( [ dx + [ dx x \ x / e\^J—x x Js x , este distribuie (numita valoarea principala a lui 1 ) 5 1 21 Exercițiu Sa se arate ca x • P- = 1 x Transformarea Fourier a distribuțiilor 99 Rezolvare Avem /x p 1 F\ = /p 1 xF\ = iim n (r-e x^(x) dx + fZ x^(x) dx\ \x ' x x iTV] = lime\0 J-, x dx ' Je x 'J = JZ, F(x) dx = oricare ar fi p £ D(R) 5 1 22 Exercițiu Fie funcția f : R* —> R, f (x) = ln |x| Aplicația f : D(R) —> C, este distribuție și = lim [ In |x| f(x) dx + / In |x| f(x) dx c1' J — , Je 1 (f) = P x ■ Rezolvare Utilizand integrarea prin parți obținem {(fY, T) = ~{f, F ) = - lime\0 (/—, ln(-x) F/(x) dx + JT ln x F/(x) dx) — lim (m^) m( £-) + f-e ^(x) dx> + fZ ^(x) — IP 1 = lime\0 - M-e) + J—, dx + Je dx) = \P x ’T> oricare ar fi f £ D(R) 5 1 23* Exercițiu Fie Th distribuția regulata corespunzatoare functiei Heaviside H : R —> R, Sa se arate ca H (x) = daca x 0 F [Th ] = -i P X+ nS Rezolvare Plecând de la relatia (Th')' = S deducem succesiv F [(Th )'] = 1 -ix F [Th ] = 1 F [Th ] = -i P1 + CS unde C este o constanta Ultima relatie este echivalenta cu (F [Th = -i + C (S, F), oricare ar fi f £ S (R) (5 1) P~, F/ x 0 1 Deoarece 100 Aplicații ale transformărilor Fourier și Laplace 2 2 g2 (F[Th], e-x2) = {Th, F[e-x2]) = {TH,^e~ V} = V fQ° e-(2) dț = 2^n jQ° e-t2dt = n §1 (2 2 \ r-e e~x e~x \ / dx + / -dx\ =0 J- x Je x J rezultă că în cazul p(x) = e-x2 relația (5 1) devine n = C (ă, e-x2) = C 5 1 24 MATHEMATICA: Definiția utilizata FM(x) = —= I' eitx Rezolvare Plecand de la relatia x • Px = 1 deducem succesiv F UPi] =2»S —iF [ix • P1] = 2n ă —i = 2n S F[px] 2''1 H + C unde C este o constanta Ultima relație este echivalenta cu F P- x = 2ni (H, p} + C(1, p}, oricare ar fi p G S(R) (5 2) Deoarece F \r1 x 1 x2 , e x x2 \ - ) = 0 în cazul p(x) = e x2 relatia (5 2) devine r 2 2 0 = 2ni / e-x dx + e-x dx JQ J-^ si conduce la C = —ni Dar 2ni H — ni si sign definesc aceeasi distributie Transformarea Fourier a distribuțiilor 101 5 1 26 MATHEMATICA: Definiția utilizata F [^](x)=■—= eitx C este continuă atunci ZT d(x) f (x) v(x) dx = {Tf, tip) T 5 3 5 Propoziție Spațiul S(R) al funcțiilor indefinit derivabile V : R —C cu proprietatea că oricare ar fi k,m G N există o constantă c^m astfel încât |xk v(m) (x) | o din S(R) converge la functia constantă 0 dacăa lim Vn = 0 n^T Transformarea Fourier a distribuțiilor 103 sup \xk T m (x) I n ' 'î' 0 oricare ar fi k,m E N 5 3 7 Definiție Spunem că o funcție liniară f : S (R) —' C este continuă daca 5 3 8 5 3 9 lim = 0 n^î lim f (spn) = 0 Definiție Prin distribuție (temperată ) se înțelege o aplicație f : S (R) —' C care este liniara si continua In loc de f(t) scriem Propoziție Spațiul distribuțiilor SZ(R) = { f : S (R) —' C | f este liniară si continuă } considerat împreună cu operațiile de adunare si înmulțire cu scalari f f + {g, t), (A = a (j, pj este un spahiu vectorial 5 3 10 Definiție Spunem despre o functie f : R —' C ca este local integrabilă daca / |f (x)| dx , oricare ar fi multimea compacta K c R Jk Utilizam notatia L1OC(R) = { f : R —' C | f este local integrabila } 5 3 11 Orice functie continua f : R Functia —' C este local integrabila dacăa 0 dacăa nu este local integrabilăa deoarece y |f (x)| dx = to f : R —' R, f (x) = J X x = 0 x = 0 104 Aplicații ale transformărilor Fourier și Laplace 5 3 12 Definiție Spunem că funcția f : R —> C este cu creștere lenta dacă există k E N și M E (0, w) astfel încăt I f (x) | ' v (, ' R : x sin x și R —> R : x cos x fiind marginite, sunt functii cu crestere lenta Functia exponentiala R —> R : x ex nu este o functie cu crestere lenta 5 3 14 Propoziție Daca funcția local integrabila f:R —>C este cu crestere lenta atunci aplicatia Tf : S(R) —> C, (Tf ,p) — f (x) w(x) dx J — este distribuție temperata Demonstrație Aplicatia Tf este liniară (Tf, av + = a{Tf, ip} + fi(Tf X) Functia f fiind cu crestere lenta, exista k E N si M E (0, w) astfel încat |f (x)| 0 este un sir din S(R) convergent la 0 atunci n^ sup |xm ^n(x)| -> 0 x€R oricare ar fi m E N si prin urmare \(Tf , 0 Transformarea Fourier a distribuțiilor 105 adică lim {Tf ,pn) = 0 5 3 15 5 3 16 Definiție Distribuțiile temperate definite de funcții local integrabile cu creștere lenta sunt numite distribuții temperate regulate sau de tip functie Celelalte distribuții sunt numite distribufii temperate singulare Propoziție (Distribuia Dirac) Aplicafia Sa : S(R) —> C, (fia, p) = p(a) este o distribuție temperata, oricare ar fi a ER Demonstrare Aplicata Sa este liniară (fia, ap + fifi) = ap(a) + fifi(a) = afia, p) + fifia, fi) Daca (pn)n>0 este un sir din S(R) convergent la 0 atunci n^^ sup |pn(x)| > 0 si prin urmare n^^ |(ăa,Pn(| = |Pn(a)| 0 5 3 17 Se poate arata ca distribuia Dirac Sa este o distribuie singulara In cazul în care a = 0 în loc de So se scrie simplu S, adica avem S : S(R) —> C, (S,p) = p(0) 5 3 18 Propoziție (Derivarea distributiilor temperate) Daca f : S (R) —> C este o distribuie temperata atunci aplicatia f> : S(R) C, (f>,p) = -{f,pfi este de asemenea o distribuție temperata, numita derivata lui f Demonstratie Aplicatia f' este liniara f ap + fifi) = -f, ap' + fifi') = -a(f, p) - fi (f, fi') = afp) + fițffi) Daca (pn)n>0 este un sir din S(R) convergent la 0 atunci sirul (pn)n>0 este un sir din S(R) convergent la 0 si prin urmare 106 Aplicații ale transformărilor Fourier și Laplace lim fpn) = - > =o- 5 3 19 Propoziție Operația de derivare a distribuțiilor temperate este liniară (af + ig)(k) = af(k) + ^g(fc) si continuă , (k) n^c (k) fn -^ f fr —> f(k}- Orice serie convergentă poate fi derivată termen cu termen oo oo £fn = f =^ = f(k) n=0 n=0 Demonstratie Dacă f, —> f atunci / „(k) \ 7 (k)\ n^C \ k I „ (k)\ / *(k) \ = (-1)k (-1)kU,¥(k)) = f (k),¥f 5 3 20 Propoziție (Multiplicarea unei distribuții cu xk) Dacă k G N și f : S (R) —> C este o distributie temperatăa atunci aplicatia xkf : S(R) C, xf,p) = (f,xM este de asemenea o distribufie temperată Demonstratie Aplicația xk f este liniara (xk f, ap + fifi) = af, xk p) + /3{f, xkfi) = a(xk f, p) + fi(xk f, fi) Daca (pn)n>0 este un sir din S(R) convergent la 0 atunci sirul (xkpn)n>0 este un sir din S(R) convergent la 0 si prin urmare lim {xkf,pn) = lim f,xkpn) = 0 5 3 21 Se poate arăta ca aplicatia df : S(R) —> C, (fif,p) = (f,dp) este o distributie daca f este distribuie si daca d apartine multimii M(R) a functiilor indefinit derivabile d : R —> R cu proprietatea ca oricare ar fi k G N exista m G N si C G (0, to) încat |d(k)(x)| C este o functie local integrabila atunci distribuia definita de functia g : R —> C, g(x) = f (ax+b) este /oo 1 yo f (ax+b) T(x) dx = — / f (x) T -o I a| J —o adicăa x — b\ , dx a x—b T f(x),T I ) / • a 5 3 24 Propoziție (Schimbarea de variabilă) Daca a,b G R sunt doua constante fixate, a = 0 si daca f : D(R) —> C este distribuție atunci funcționala f (ax+b) : D(R) —> C, j (ax+b),T(x)) = I1 f (x),T x—b a este de asemenea distribuie 5 3 25 Cazuri particulare importante: Translatia : Omotetia : j(x+b),T(x)) = (/(x),T (x —b)) 1x (f (ax), T(x)) = — ( f (x), T ( - ) |a| a 5 3 26 Exercițiu Sa se arate ca b(x — b) = 6b si b(ax) = pî-Ț 5 |a| 5 3 27 Definiție Spunem despre o distribuție f gD'(R) ca este nula pe o mufiime deschisă D c R daca (f, T) =0 oricare ar fi t cu supp t C D Complementara celei mai mari mulțtimi deschise pe care f este nulăa se numeste suportul lui f țsi se noteazăa cu supp f Spunem căa f este distributie cu suport compact dacăa supp f este mulțtime maărginităa 108 Aplicații ale transformărilor Fourier și Laplace 5 3 28 Distribuția 6a este cu suport compact (punctual) deoarece supp6a = {a} 5 3 29 O functie f : R —> C absolut integrabilă, adică astfel încât Zo |f (x)| dx C este o distribuție temperata atunci aplicația F[f] : S(R) C, (F[f],p) = (f,F[^]> este de asemenea o distribuție temperata (transformata Fourier a lui f) Demonstrare A se vedea 5 3 31 Propoziție Transformarea Fourier a distributiilor F : s'(r) -^ s'(r) : f F[f], (F[f],^ = (f,F[^]> este bijectiva și continua Inversa ei este F—1 : S'(R) -^ S'(R) : f F—1[f], (F—1[f = (f, F—1 [^]> Demonstratie Avem si (F[F— 1 [f]],^( = (F— 1[f], F[^]> = (f, F— 1[F[^]]^> = (f,V) Transformarea Fourier a distribuțiilor 109 = (F[f], F-1 f F[F-1 = {/,?) oricare ar fi G S(R) și f G SZ(R) Dacă fn —> f, adică dacă {fn, f, fi) oricare ar fi G S(R) atunci = {fn, F[^]> {f, F[^]> = {F[f],fi) și prin urmare fn > f F[fn] > F[f]- 5 3 32 Din relația (5 3) rezultă ca transformarea Fourier a distribuțiilor poate fi privita ca o prelungire a transformarii Fourier a functiilor 5 3 33 Propoziție Relațiile (F [f])(k) = F [(ix)k f ] F [f(k)] = (-i£)k F [f ] au loc oricare ar fi distribuția temperata f Demonstrație Utilizand proprietătile transformarii Fourier a functiilor obtinem {(F[f])(k),fi) = (-1)k {F[f],fik)) = (-1)k f F[^fc)]) = (-1)k f HO" F[^]> = = (F[(i£)k f],fi) {F[f(k')],fi) = {f(k), F[^]> = (-1)k f (FM)F) = (-1)k f F[(ix)k = (-1)k {F[f], (ix)k fi) = {(-ix)kF[f],fi) oricare ar fi G S (R) si f G SZ(R) 5 3 34 Transformarea Fourier joaca un rol fundamental în matematica si aplicabile ei Integrala / ei?x xk dx J —^ nefiind convergenta pentru k G N, rezulta ca funcțiile constante si cele poli-nomiale nu admit transformate Fourier dacăa ne limităam la abordarea clasicaă Utilizarea distribuțtiilor temperate largețste considerabil posibilităațtile de folosire a transformarii Fourier în diverse modele matematice 110 Aplicații ale transformărilor Fourier și Laplace 5 3 35 Propoziție Daca șirul de funcții local integrabile (fn)n>0 converge uniform pe fiecare mulțime compacta la funcția f atunci și prin urmare lim Tfn = Tf n^c J J oricare ar fi k = 0,1,2, n^c Demonstrație Pentru orice GD(R) exista r> 0 astfel încat supp C [—r, r] Deoarece w este funcție marginita avem fnw “———]> fw si prin urmare pr fi r rr fn(x) w(x) dx = / f (x) w(x) dx = {Tf ,ipf rn r f r 5 3 36* Exercițiu Seriile de distribuții oo V einx si C : w ( Y'' §na, \ = lim (5^ Sna, = lim ip(na) n=0 n=0 n=0 Transformarea Fourier a distribuțiilor 111 este bine-definită și liniară Dacă —> 0 în D(R) atunci există r > 0 astfel încât suppC [—r, r], oricare ar fi m = 0,1,2, I (k) / J m^c sup \ 0, oricare ar fi k = 0,1,2, x€[-r,r] Daca nr G N este astfel încat nra m(na)| 0 Convergența celei de a doua serii rezultă din relația oc oc oc y? ^na ^na + ^n(-a) ■ n=-c n=0 n=1 Figura 5 6: Graficul funcției f 5 3 38* Teoremă In spațiul distribuțiilor D'(R) are loc egalitatea oc R functia periodica cu perioada 2n definită prin relația 2 X X f (x) = 2 — — pentru 0 C | f este liniara si continua } Tf : d(r2) —> C, {Tf C, T) = T(a b) f : D(R2) C, (f =- {f ) : d(r2) ► C {^y= {f, G^y , etc- 5 3 42 Propoziție Daca f g:D(R) —>C sunt doua distribuții atunci funcționala f xg : d(r2) —>C f xg,T) = f (x) {g(y),T(x,y)}] este o distribuție, numita produsul direct al distribuțiilor f si 9- 5 3 43 Din relatia {6a x âb T) = {6a(x), {6b(y),T(x,y)')') = {6a(x),ip(x,b)) = p(a,b) rezulta ca 6a x 6b = 6(a,b) ■ 5 3 44 Definiție Fie f g: R —> C două functii astfel încât integrala Zo f (t) g(x - t) dt este convergenta oricare ar fi x G R Funcția Z f (t) g(x - t) dt se numeste convolutia functiilor f si g 5 3 45 Din relatia ZOO f (t) g(x - t) dt = f (x - u) g(u) du = (g*f )(x) -O J — O rezulta ca f *g = g */■ 5 3 46 Definiție Spunem ca sirul (yn)n>o din D(R2) converge la 1 în R2 daca: a) oricare ar fi r> 0 exista nr G N astfel încât pentru n > nr avem dn(x, y) = 1 dacă x2 + y2 (x + y)} n^^ exista pentru orice functie GD(R) si nu depinde de sirul (nn)n>0 convergent la 1 în R2 ales atunci functionala f *g : D(R) —> C, f *g, = nljmo(f (x) xg(y), pn(x, y) o c D(R2) converge la 1 în R2 atunci lim nn(x, b) C se poate prelungi unic pana la o distribuie g : S (R) —> C utilizând relația {g,?} = unde n GD(R) este o funcție aleasâ astfel încât Transformarea Fourier a distribuțiilor 115 n(x) = 1 pe o mulțime deschisă ce include supp g Se poate arătă că F[g] este o distribuția temperată regulată asociată funcției F [g]« ) = (g(x),n(x)ei«x) si apartine spatiului M(R) (a se vedea pag 106-21) 5 3 50 Teoremă (Transformarea Fourier a convolutei) Daca f ESz(R) este o distribuție temperată și daca g gD'(R) este o distribuție cu suport compact atunci F[f * g]= F[g] -F[f ] Demonstrapie A se vedea 5 3 51 Propoziție Daca f : R —> C este local integrabila cu crestere lentă §i Fr : R —> C, fr (x) = f (x) 0 pentru pentru |x| r atunci a) b) c) lini, x Tfr = Tf limr^rc> F[Tfr] = F[Tf] F[f ](£)=«& Jei«x f (x) dx în S '(R) —r r Demonstrație a) Oricare ar fi G S(R) are loc relatia Zr r oc f (x) w(x) dx = / f (x) w(x) dx = (Tf ,F) r J- ]■ c) Relatia de la b) se poate scrie succesiv sub formele 116 Aplicații ale transformărilor Fourier și Laplace limr^tt F[Tfr ], = F[Tf], limr^tt (Tfr ,FM) = (Tf,FM) r tt lim f f (x) F[^](x) dx = f f (x) F[^](x) dx r^tt—r —tt r tt tt tt ^m f f (x) f e1?M£) df dx = f f (x) f e1?x^«) df dx tt —tt lim r^tt r —tt r lim f e1?x f (x) dx ^(f) df = f f e1?x f (x) dx ^(f) df r^tt —r r fe1?x f (x) dx, ^(f n = (f [f ](f), • —r tt tt — tt tt —tt —tt 5 3 52 Dacă nu distingem doua functii care difera doar pe o mulțime de masura nula, spațiul funcțiilor de pătrat integrabil (Lebesgue) tt A y |f (x)|2 dx —tt considerat împreună cu produsul scalar L2(R) = C tt țf,g) = j f (x) g(x) dx —tt este un spațiu Hilbert 5 3 53 Teoremă (Plancherel) Transformarea Fourier F : L2(R) —> L2(R) definită prin r F[f](f) = lim y e1?x f (x) dx în L2(R) —r este bijectivă, bicontinua si are loc egalitatea Parseval [F[f], F[g]) =2n (J,g) oricare ar fi f,g G L2(R) 5 3 54 In cazul particular f=g, egalitatea Parceval devine ||F[f]||2 = 2n ||f ||2, adică tt tt y |f [f]o2 =2n /f (x)|2 dx- —tt —tt 5 3 55 Teoremă Daca f,g G L2(R) atunci în S'(R) are loc egalitatea F[f * g] = F[f ] -F[g] Demonstrație Dacă f, g G L2(R) atunci F[f], F[g] G L2(R) Din relațiile [J |F [f ](C)F [g]O< < J |F [f J |F [g](C)|2 < [ff |f (t)g(x-t) ^(x)| dtdx\2 < [ff |f (t)|2 c c dtdx] [ff |g(x-t)|2 |^(x)| dtdx] <||f||2 ||g||2 [JJ |^(x)| dtdx]2 < rezultă că F[f]-F[g] și f *g definesc distribuții temperate Oricare ar fi gS(R) (F[f*g],^) = f *g,d7[^]) F[^](xHf(t)g(x-1)dtdx = f f (t) f g(x - t) F[^] (x) dx dt = f f (t) f F[g(x - t)] (O d< dt = f f (t) f F[g] (e) ^(C) ei?td£ dt = f F[f]«) F[f]«) <) d< 117 118 Aplicații ale transformărilor Fourier și Laplace Capitolul 6 Transformarea Laplace a functiilor 6 1 Exemple 6 1 1 6 2 Aplicații 6 2 1 6 3 Definiții si rezultate fundamentale 6 3 1 119 120 Aplicații ale transformărilor Fourier și Laplace Capitolul 7 Transformarea Laplace a distribuțiilor 7 1 Exemple 7 1 1 7 2 Aplicații 7 2 1 7 3 Definiții si rezultate fundamentale 7 3 1 121 122 Aplicații ale transformărilor Fourier și Laplace Capitolul 8 Transformarea Laplace finită 8 1 Exemple 8 1 1 8 2 Aplicații 8 2 1 8 3 Definiții si rezultate fundamentale 8 3 1 123 124 Aplicații ale transformărilor Fourier și Laplace Bibliografie R J Beerends et al , Fourier and Laplace Transforms, Cambridge University Press, 2003 R N Bracewell, The Fourier Transform and Its Applications, Mc Graw Hill, 2000 D J Brady, Optical Imaging and Spectroscopy, Wiley, 2009 N Cotfas, Elemente de Algebra Liniara, Editura Universității din Bucuresti, 2009 N Cotfas si L A Cotfas, Elemente de Analiza Matematica, Editura Universități din Bucuresti, 2010 M Fox, Quantum Optics An Introduction, Oxford University Press, 2006 J F James, A Students Guide to Fourier Transforms, Cambridge University Press, 2011 J Kauppinen, J Partanen, Fourier Transforms in Spectroscopy, Wiley-VCH, 2001 V S Vladimirov, Ecuațiile Fizicii Matematice, Editura Științifica si Enciclopedică, Bucuresti, 1980 V S Vladimirov si alți, Culegere de Probleme de Ecuațile Fizicii Matematice, Editura Stiintifcăa si Enciclopedicăa , Bucuresti, 1981 http://fpcm5 fizica unibuc ro/ ncotfas/ 125